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ABSTRACT 


The theory of excitation and charge exchange in proton-hydrogen 
collisions is discussed within the framework of the impact parameter 
method. Consideration is given to the importance of proper boundary 
conditions. The time dependent equations linking the amplitudes of 
the target with those of the rearranged system are written in matrix 
form. By eliminating the rearrangement amplitudes from these equa- 
tions a second order matrix equation is derived which may he used as 
a basis for successive approximations which are automatically second 
order. The theory is generalised with the aid of a two center 
expansion of the electronic wave function. The method is illustrated 
by computing first and second order distortion approximations for the 
reactions H + + H(ls) - H + + H(2s) and H + + H(ls) - H)ls) + E + and the 
results are compared with previous calculations. 


51 Introduction 


In this paper we shall he concerned with methods of calcula- 
ting the e xcit atio n and charge exchange cross sections of proton- 
hydrogen collisions. Our interest will lie mainly in the kilovolt 
region and for this purpose it is sufficient to use the well known 
"impact parameter method"' 1 ' (IPM) in which the protons are 
treated as cl as sical - , particles mov ing with constant relative 
velocity. Thu s only the electron need he treated quantum 
mec h a ni c all y as it moves in the time-dependent field of the "in- 
finitely massive" protons. 

Itofoirtunately, the apparent simplicity of this problem is 

plagued hy the usual bugbear of atomic scattering theory; the 

necessity of accounting for the infinite number of hydrogenic 

states. Specifically, the strong coupling, found in calculations 

based on expansions in atonic orbitals seems to indicate a strong 

coupling with states excluded from the calculation and the sub- 

(o) 

sequent slow convergence of this approach . This difficulty 

is somewhat enhanced by the fact that virt ually all inelastic 

channels are energetically permissible. However, relatively 

little is known about this coupling and de taile d investig a ti ons 

(3) 

are only now being made . At high energies recent studies 
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indicate that the Influence of the continuum may veil he decisive 
Moreover, since it is not clear hov this latter effect could he 
incorporated in the usual expansion methods, it is possible that 
these methods may converge to an incorrect result. In order to 
investigate this possibility ve have devised an appr oxima tion 
scheme based on an expansion in atomic orbitals hut vhich simul- 
taneously includes an effect firom all those states (including 
the con tinuum ) vhich have not heen explicitly included in the 
calculation, ihus, to parellel the usual "first order methods" 
ve derive a set of "second order methods". 

A partial review of first order methods is given in Section IH 

using a matrix notation introduced in Section U. In contrast 

(5) 

with previous authors ve place particular emphasis on the 
importance of correct boundary conditions. For resonant charge 
transfer ve introduce a distortion approximation and a modifi- 
cation of the Br i hkman -Kramers approximation both of which have 
the correct bo undar y conditions. However, the main purpose 
of the review is to show clearly the analogies between the first 

\ 

'order methods and the second order methods introduced in Section 
IV and to present the matrix equations necessary in the derivation 
of these second order methods. Finally, in Section V, we calculate 
cross sections for the reactions 


H + + H(lS) H(1S) + H* 


f 


and 


H + + H(1S) - H + + H(2S) 


according to the distortion approximation in both the first 
and second order methods and discuss the significance of the 


results. 


$11 Notation 


Let e, A and B de n ote the electron, the target proton and 
the Inci dent proton, respectively. Let R he the position vector 
of B r elat ive to A and let r., r_. and r he the position vectors 

“A “33 — 

of e relative to A, B and the mid point of AB- In the IPM it 
is assumed that A remains fixed while B moves in a straight line 
•with a constant speed, v. Thus R = p. + vt where t is the time, 
chosen such that at t = o the protons, A and B, have a ” -imum 
separation, j j , which is the impact parameter for the 
collision. 

(6) 

The time-dependent Schroedinger equation (in atomic units } v 
for the complete electronic wave function, Y (r, t), is 

(2. + V B ) V = o (la) 

or 


<*3 + V V - ° 


(lb) 


where 


2 


A = 




at 



(2a) 




i V 2 + 




B 


(2b) 


1 » ? 

I 

• 
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» 

1 

1 

1 

If ^ (r) denotes a hydrogenic eigenfunction with eigenenergy 

£ so that 
n 



(i V 2 + — + e ) 0 (r) = o 

( 5 ) 


r n y n - 

then 


1 

1 

1 

[ 

i 

§ A n 6X33 ^ _i ^ - - + ir v®t + e n t} 

( 4 a) 


=0^(23) exp c+i {§!*£. - ir v2fc 

(4b) 


are exact solutions of the unperturbed equations 

T a =0 

( 5 a) 


A n 


* 

0 

11 

(5b) 


and form two mutually exclusive complete orthogonal sets. 



i 
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Consequently we may expand Y as 


Y = A 


(6a) 


or 


Y 




(6b) 


where 

A 


and 


B 


are row matrices with elements 


B 


$ _ and 5 while A and B are column matrices with 
n n 

elements a^(t) and b^Ofc.) • The scalar products (6a, h) imply 
integration over the continuum as well as summation over all 
hound states. 


if (la, h) is solved subject to the boundary condition 


( 4 ) 


t -» - “ 


£ (vB - V*tf l/V 5 


np 


(v) 


then the probability of excitation of the target from an ini t ial 
state (p) to a final state (q) is 


(P, q) = Lira \ a. Q (i) | 2 
t - + ® 

and the probability of charge transfer to a state (q) is 


(8a) 


? B (P, <l) = Lim jb (t)j 2 

~ t -• + to ^ 


(8b) 


The corresponding cross sections are obtained by integrating 
over all possible impact parameters 


\ B (P ’ a) 


d £_ 


'A, B 


(p, q) 


§ IU First Order Methods 


A 3 

■By projecting the states § and § on (la) and (lb) ve 


obtain 


(§ a |t a Iy) - - (s A jv B j>) 


($ B J^|y) = - (sVjy) 


vhere 


(^ B n U 1*) = ? A,B n (r,t)x (r,t) Y(r>t) 


_A,B 


Equations (10a, b) are entirely equivalent to (la,b) and form 
a convenient starting point for our discussion of first order 
methods . 

Substituting (6a) in (10a) and (6b) in (lOb) gives 

i -i uXi^A 


(10a) 


(10b) 


( 11 ) 


(22a) 


§ -i (s b |v a U b )b 


(12b) 


* %■ 
» 
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A A B 

where (§ IVhls ) and (§ |V, js ) are square matrices whose 
(n, m)^ elements are ( Jv. j § A ) and ( § 3 jV ) ). Approxi- 



nations to (12a) may he made by retaining only a few specific 

1 


elements of (§^jYgi^)* for example, by retaining only the 

i 


diagonal elements the equation is uncoupled and solved by 

i 

- 

a n = 6 i® exp { i I < 5& J V Bl‘ A n )ar '} 

f 

(13) r 

‘ 

! 

which satisfies (7) • To obtain an estimate of (q £ p) we 

neglect all terms containing states other than (p) and (q). 

» 


to give 

t 

I 

| 

a = V a + V a 

( 1 ^) ; 

i 


q qq 0 qp p 

| 

' 

\ where 


i - 
1 ' 

A 

i 

\ • 

v = (? a J v b U b J 

nm n ' ±s 1 m 

(15) t 
♦ 

I 

\ 

m 



Equation (l4) is solved. •by^T) 


a = exp(i f V dt‘"V f dt* V a exp -f-if V dt"l 
4 l J 44 J ~_ a qp p 1 J 44 J 

which, from ( 8 a), leads to the excitation- probability 

, CO ^2. 

V 1 ' a) 51 1 J./* V p 6X11 I' 1 J V 44 “'} I 

Substituting (U) for a p in ( 17 ) gives the well known distortion 

(7} 

approximation introduced by Bates w ' in his calculations of 
Q A (ls, 2s) and Q^(ls, 2p). Approximation ( 17 ) will prove useful 
when further methods of computing a p are discussed later in this 
paper. 

To procede further with (12a) we may solve numerically 
the coupled equations which result when certain specific ele men ts 
are retained. For example, Q (is, 2s) has been calculated in the 
1 sA/ 2 sA^^ and in the lsA/ 2 sA/ 2 pA^^ approximations to (l 2 a). 
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By including higher states more accurate results may he ex- 
pected. However, the convergence is likely to he slow since no 
allowance has heen made for coupling with rearrangement states. 

Moreover, since the boundary condition (7) is embedded in the 
continuum of B the usefulness of (l2a) is greatly restricted^ ^ ; 
indeed from a computational point of view it is virtually useless. 

As an alternative to ( 12 s, h) we interchange the substitutions 
of (6) in the right hand side of (10) to obtain 

dA A b 

= i ({ |V 3 U ) B (lBa) 

fr =i U E :V A |t A ) A (18b) 

Clearly, if we already have a knowledge of A, (l8b) may he 
integrated to give the charge exchange amplitude 

CO 

p b (p, a) » | f at (§ B !v a |i a ) a ] 2 (19) 

—CO * 

She well known Brinkman -Kramers ^ ^ (EK) approximation is ob- 
tained from (19) by making the substitution a = 6 . Ihis is 

n np 
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inconsistent with ( 7 ). Instead of 5 we substitute 

np 

6 (vR - v^t) 17 ^ for a in ( 19 ) to obtain a modified Br inkman - 
up n 

Kramers (MBK) approximation which is consistent with ( 7 ). She 

situation may be further improved by substituting ( 13 ) for a 

n 

in ( 19 ) to give a distortion approximation analogous to that des- 
cribed above for excitation. 

As a basis for successive approximations, however, (l3a, b) 
is defective. This is illustrated by considering the diagonal 
approximation 


a = i ( § A |V_ j f ) b 
n n 1 £ ! n n 


(20a) 


= i U B iv. !$ A ) a 
n 1 A' n ; 


(20b) 


which has the general solution 


& n = c i e 1Y n + c 2 


b^ = c x e lY n - c 2 e lY n 


(21a) 


(21b) 


where 



and ci and c 2 are arbitrary constants. It is easily shown that 
Y n vanishes rapidly for large negative t so that (21a) is in- 
consistent with ( 7 ) • This fundamental discrepancy cannot be 
resolved by the inclusion of further bound states. 

To overcome the defects of (12a, b) we .expand Y in an 
over complete set 


■* « ^ A + $ B B 


where A and B are to be determined by (10a, b). The "beauty" 
of this two center expansion is that it makes explicit allowance 
for each reaction path and thus circumvents the defects of the 
single center expansions (6a, b) where rearrangement states are 
awkwardly contained in the continuum and thus confused with 
ionization states. Perhaps an even better expansion would re- 
sult if, for example, the second term on the right hand side of 
(23) were restricted to a summation over bound states only. We 
would then have a clearer physical interpretation for the continuum 
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elements of A which would correspond to pure ionization. This 
introduces a slight compli cation, however, and since we are here 
mainly concerned with excitation and charge t rans fer we shall 
simply exploit the symmetry of (25). Clearly the hound state 
coefficients of (25) mast co-incide with those of (6a, h) at 
infinite proton separation and there is no need to alter the pro- 
bability definitions (8a, b) 

Substituting (25) in (10a, b) gives 

5T * U A |5 B >f|- iU A |v B |s A ) A + l(» A |V A |t B ) B 



(i B | A t£- *J* B |V A |* B ) B + lUVsl* 4 ) a 


Again it is instructive to consider the diagonal approximation. 
This approximation co-in cides with the two state approxxmation 
of Bates ^ ^ and has been solved by McCarroll^^^. 'Pho solution 
may Be put in the form 

a n " 6 np {vS ‘ v2t) ” i ^ V e 3 n 

r 

- I ’ 

b n “ 6 n^> * vR “ e±a »sin 


(24a) 


(24b) 


( 25 a) 


(25b) 


where a and B are real bhA vanish as t •bends to - ®. She bound- 
ary condition (7) is therefore satisfied by (25a) ia contrast 
with (21a). , 


S IV Second Order Methods 


Broadly speaking, the methods discussed so far may be 

considered in the following way. By restricting the numb er of 

available states to a select few we obtain a tractable model 

which may be treated precisely by numerical methods. tmi* is 

equivalent to allowing for an infinite number of transitions 

between a limi ted number of states and it neglects completely 

effects due to all ot h e r inelastic processes . Thus the continuum, 

which may play an extremely important role in the intermediate 

stages of the collision, is neglected. Methods which attempt 

to take account of such transitions axe loosely termed "second 

order". Far example, the second Bora^ 1 ^, the impulse and 

( 4 ) 

the continuum distorted wave approximations may be thought of 

in t his l i g ht * The significant feature of these approximations 

is that, far resonant charge transfer form the (is) state, they 

( 15 ) 

all predict the high energy behaviour 


Qgds, 1 b) ( 0 . 29 ^ + 'u v) (is. Is) 


(2 6) 


where p is a small (and for our purposes, insignificant) constant, 
and Q^-g(ls, Is) is the Brinkman-Kramefs cross section 


1 


( 27 ) 


17 


Q^gCls, Is) o (1 + iv 2 ) (-nao 2 ) 

Thi s is in very sharp contract with the high energy behaviour 
of first order methods. Hie two state Bate s -Me Carroll approxLxna- 

W f V) \ 

tion behaves Hire 0 (is. Is) ' at high energies and there is 
no reason to suspect that this will be significantly altered by 
the inclusion of further bound states. However, the second order 
methods mentioned above give poor results at low energies (below 
50 kev) and seem to he incapable of systematic improvement . More- 
over one's confidence in (26) is, to some extent, undermined by 
the likelihood that the Born series diverges for rearrangement 

(16) 

collisions' . Clearly it would be advantageous to devise an 
approximation scheme which includes effects from all inelastic 
processes and which may simultaneously ‘be treated by suc^ssive 
approximations based on the inclusion of the more significant 
discrete states. We now proceed towards this end. 

First, let us consider how we could, in principle, solve the 
coupled equations (l8a, b) . Suppose we start by substituting 
a given first order approximation to A in (lBb) . He could then 
calculate each element of B and substitute the result in (l8a) . 
*tm « would involve an infinite number of coefficients, b^. 


I 


corresponding to the infinite number of available hydrogen! c 
states. Moreover, the most important of these (corresponding 
to rearrangements) are embedded in the continuum. Nevertheless, 
having surmounted these barriers ve would now perform the in- 
finite sum implied in the R.H.S. of (l8a) ana proceed to 
calculate the second order approximation to A." We would then 
repeat this process and hope that A would converge in successive 
Iterations. However, if we could somehow eliminate B from 
(lBa, b) the intermedi a te stages described above would automatically 
be contained in the resulting equation for A. This is the cen- 
tral idea of the present paper. The process of eHminft - t-.-Sng b 
from (l8a, b) is now performed. 

Assuming the existence of the matrix ($ B |V B |§ A ) we consider 
the product 


($ B jv“3[i A )(i A lv B |« B ) - (s B |v"^v B j$ B ) - - I 


where X denotes the unit matrix. In the derivation of (28) 

A 

we have made use of the closure property of } together with 
the orthogonality of $ . Similarly, we can show that 



-is - 


(» B |t‘ 1 b |I A ) -2 


(29) 


bo that, formally. 


< iv _i B i - uv^r 1 


(50) 


This enables us to vrite (18a) as 


B-.-i (^Iv'^j# 4 ) A 


(51) 


and, differentiating both sides. 




( 52 ) 


Making use of (5b) and assuming the validity of Green's Theorem, 


ve have 



(33) 


1 ft <«V 1 :b I« A > - 

whence (32) becomes 



<*V l B |**)g§- 1 (i?l^-i B | S A )f£ 


(34) 


^■n 

We now substitute (l8b) forr^in (34) and • pre -multiply by 
A B 

( ? jV B j$ ), making use of (30) and the closure property 
B 

$ , to obtain 



U A |V B I B V" 1 B lt A fr+ (t A |V B V A | 4 A ) A . o 


(35a) 


\ which is the desired result. Thus, ve have replaced two coupled 
first order matrix equations by a single second order matrix 
equation. A similar analysis, in which A is eliminated from 
(iBa, b) yields 


d^B 

dt 2 


- 1 <* B IWAI^§ + (sVa^I* 8 ) b - ° 


(35b) 
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If we now perform calculations in which only a limited 

A 

number of elements of § are retained in the matrices 
<$ A |V B T B V” 1 iJ | 5^ and ( * A |V b V a | $ A ) we already include continuum 
effects due to the intermediate transitions of (l8). Eq. (55a) 
is most useful for calculating excitation amplitudes while (55b) 
suffers from the difficulty of including the proper boundary 
conditions and in this respect it resembles (12b) . The usefulness 
of (35 a * b) is therefore restricted. To overcome this limitation ' 
we must use (24a, b) as our starting equations rather than (l8a, b) 
In this case the derivation of the second order equations is 
somewhat more complicated although the principles involved are 
similar to those employed above. A simpler derivation, which 
has the additional advantage of avoiding assumptions about the 
existence of the inverse matrix, follows from the identities 

<* A |W" 1 B T Al' r) = ( * A |V A l’ f) 


which may easily be verified with the aid of (la, b). Equations 
(36 a, b) may be considered as the starting point of the second 
order theory replacing (10a, b) of the first order theory. For 



(36a) 


(56b) 


example, by substituting ( 6 a) in ( 36 a) and ( 6 b) in ( 36 b) ve ob- 
tain (35a) and (35b), which are the second order analogs of ( 12 a) 
and (12b). Similarly by substituting (23) for f in ( 36 a, b) 
ve obtain the second order analog of (24a, b). The result is 

0 - * (» A IW‘ sl* A ) -uV/jji/) a 

- - 1 ff 

+ (* a IW'Va - VaI* b 


- i (^IVAV'hlS 8 )! * (i B |v A v B h B ) B 


(37a) 


- - «V0* i ^ 


(37b) 
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Finally, we note the following useful relationships which may 
easily be verified 


uVbVs'i *) (j*|^- V B|x) + (* A |iv| + (i+Uv b )t a |x) 

(t °l V A T B V rU) - -lit +U B |4V^ + (l+tV A )T B |i) 



1 




V. 


(38a) 

(38b) 


§V Detailed Calculations 


To illustrate the method we calculate a (t) by applying 

is 

the diagonal approximation to (35a). The solution may he re- 
garded as the second order analog of ( 13 ) and is used to cal- 
culate inelastic cross sections by substituting it into the 
appropriate first order matrix elements. This procedure is 
slightly inconsistent and the results obtained are not truely 
representative of the second order approach. For example, by 
applying the diagonal approximation to ( 37 a, b) we would have 
obtained a second order equivalent of the two state Bates -McCarr oil 
approximation. However, such a calculation would be rather complex 
and will therefore be delayed to a later date. 

Our equation for a is (t) is 


^5- a xs (t) + (P ■ i Q) axs (t) + H a i s (t) = o 


V 

\ 


where 



- 25 “ ~ 

? “ " ft ^ A is| ?^, v b I$ A is) 

a - ($ A xs|v b -v a + Jho) 


2 - ($ ti v A v BiV 


It toast be solved subject to 


aig(t) - (vR - v 2 t)“ i ^ v 
*fc - 00 

X 


( 41 ) 


Details axe given in the appendix. We denote the result of this 
second order calculation by & lB ( t) to distin g uish it from the 

corresponding first order result, (U), -which we denote by a^ ls (t). 

Figure 1 shows ' Jaxs\t)j for p = 0.13Q5 and v=CL. Note 
the deep trough which occurs just after impact and the slight 
decaying oscillation which follows as the interproton distance 


( 2 \ 

increases and ja v isC^I settles to a constant, somewhat 

less than unity. This behaviour is typical, the depth of the 

trough decreasing with incleasing impact parameter and velocity. 

|a^ ls (t)| is, of course, every where unity since, unlik e 

(a) 

]a' y ig(t)|, it does not share probability with competing in- 
elastic channels. 

First and second order distortion approximations to the 
excitation probability, P^(ls, 2s), are obtained by replacing 
ais^) ( 17 ) °y a^ ls and a^ 2 | s and the corresponding cross 
sections, ft Dl A (ks, 2s) and Q^^Cls, 2s), are tabulated, together 
with the first Bom cross section, in table 1. Q^ 1 (is, 2s) is 

xx 

smaller than Q A (ls, 2s), especially at low energies. This re- 
sults from, the presence of the phase factor in the overlap integral 
of the former. Q^^Cls, 2s) is smaller than Q Dl A (ls, 2s) as a 
consequence of the allowance for inelastic processes inherent 
in a x s (t). fhe comparative purposes, the above results are 
shown in figure 2 together with some calculations of previous 
authors. It is seen that the distortion approximations are in 
broad agreement with the 1sA/2sA approximation to (l2a) and the 
lsA/2sA/lsB approximation to (24a, b) but show a marked dis- 
agreement with treatments which take account of excitation to 
the 2p levels. ' *' where the computed cross sections are 
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B 

larger than Q ^(ls, 2s). This is not too surprising, for the 
optically allowed ls-2p transitions are expected to have large 
cross sections, and since the 2p and 2s states are degenerate 
we may expect a strong coupling "between these states and a conse- 
quent I n c rease in the computed cross section, Q^(ls, 2s). How- 
ever, as ve have seen, the second order treatment has the effect 
of reduci n g the cross section in the distortion approximation 
and if a similar effect should occur in a second order cal c ulat ion 
With allowance for intermediate 2p states it appears likely 
that the fi na l result may "be close to that of the first Born ap- 
proximation. 

For resonant charge transfer ve calculate the MBK and ais- 

MBK T) 

tortion cross sections, Q ^(ls. Is) and Q ^(ls. Is), described 
in Section 111. The second order distortion approximation is oh- 
tained by substituting 6 , a' ' (t) for A in (19). The resultant 

H^XS XS 

cross section, Q^ 2 (is. Is), is tabulated together with (is, Is), 

■D B 

BK D 

Q ^(ls, Is), Q ^(ls,!^) and the Bates -McCarroll cross section 
BM 

Q B (ls, is) in Table 2. Of the above cross sections only 

•RFC 

H us. Is) is derived with incorrect boundary conditions for 
aig(t) and, as may be seen from the table, it is in considerable 
disagreement with the other tabulated cross sections over most of 
the energy region considered. The difference between (is, Is) 

and Q^gCls, Is) illustrates the distinction between first and 

0^2 (is. Is) is in extremely good agreement 


second order effects. 



BM , . 

with Q ^ &s. Is) and this seems to indicate that the hack 
coup ling effects of Q ^(ls, Is) are largely accounted for in 
Q^tls, Is). However, it should not he concluded from this 
that these cross sections provide an accurate representation of 
the exact result since the methods considered have taken no 
account of contributions from terms other than a ls (t) . 

Such contributions, if important, should he appar ent in a 
proper two state calculation based on (37a, b). 
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Appendix 

Ehe coefficients of (40) are found to "be 

Qsl + e^-jd-e Kl (A2) 

H - 5 (1 - e' 23 ) (A3) 

•with 

Ko = e 23 E(2R) - e" 23 E (-2R) (A4) 

K x . |Ko - R [e^CaR) + e _2R S (-2R)] (A5) 


I 
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E(x) = J ^ e~ 7 dy 


Writing 


&1S “ a e 


ve have 


£ + pa + (H - a Qa) a = o 


a + (j 4 + p) a » P | 


vhere 




a = 9 


(A1G) 


The boundary condition ( 7 ) implies that as t — - » a -* 1 
and a -» g - • With these considerations in mind it was found con- 
venient to introduce 


y = 


a 2 H 


(AH) 


where 


X = ^ (Ko + Kx) 


(A12) 


and to replace (a 8) and (A 9 ) by three coupled first order 
equations 



(A13) 


33 


£ 


- RaE e X (A14) 


(A15) 


with the "boundary conditions 


(A16) 
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Cross sections for H + 

+ H(ls) - E + 

+ H(2s) in units 

Of TSo 2 

leg (Energy (keV) ) 

Q^(is,2s ) 

Q^(ls,2s) 

Q? 2 (ls,2s] 

1.00 

5 .37-10 

6. 79x10 ~ 2 

2.5Sxl0“ 2 

1.25 

4-. 10x10“ 1 

1.31X10 -1 

5.70xl0“ 2 

1.50 

2.7610“- 

l.i-oxlO -1 

B.OExlO * 2 

1-75 

1.72x10*- 

1.21x10 ~ 

7.0 ^-' 0* 2 

2.00 

1.05x10 _L 

8.42x10 2 

6. 07x10 ~ 2 

2.50 

3. 42x10 “ 2 

3- 22x10 “ 2 

2.57xl0“ 2 


Table 2 


Cross sections for H + + H(ls) -* E(ls) + H + in -units of nao 2 


Energy (kev) 

5 

25 

50 

100 

400 

1000 

^(ls,ls) 

2 . 01 X 10 2 

l.SSxlO 1 

3.37 

4 . 0 CX 10 -1 

1.02xl0~ 3 

7.94x10"' 

0 ^( 13 , 18 ) 

5-02 

3.81 

1.19 

1 . 83 X 10 ' 1 

6.15x10"* 

6.87x10“ 

q£ x (1s,1s) 

4.38 

4.87 

1.52 

2.3LS10' 1 

7.86x10"* 

8.77x10"' 

0^ 2 (ls,ls) 

4.08 

3.08 

1.00 

1.55X10" 1 

5.29x10"* 

5.5xio" 6 

Q^ M (1s,1s) 

11.3 

3.13 

7.70x10* 

’ x 1.15x10"- 

5.38x10"* 

5.82x10"' 
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